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Abstract. We present high-dispersion spectra of two extremely massive star clusters in galactic merger remnants, 
obtained using the UVES spectrograph mounted on the ESO Very Large Telescope. One cluster, W30, is located 
in the ~ 500 Myr old merger remnant NGC 7252 and has a velocity dispersion and effective radius of cr = 
27.5 ± 2.5 km s~^ and RcH ~ 9.3 ± 1.7 pc, respectively. The other cluster, G114, located in the ~ 3 Gyr old 
merger remnant NGC 1316, is much more compact, Rga — 4.08 ± 0.55 pc, and has a velocity dispersion of 
a — 42.1 ± 2.8 km s^^. These measurements allow an estimate of the virial mass of the two clusters, yielding 
MdyniWSO) = 1.59(±0.26) x lO'^A^e and Mdyr^iGlU) = 1.64(±0.13) x lO'^X©. Both clusters are extremely 
massive, being more than three times heavier than the most massive globular clusters in the Galaxy. For both 
clusters we measure light-to-mass ratios, which when compared to simple stellar population (SSP) models of the 
appropriate age, are consistent with a Kroupa-type stellar mass function. Using measurements from the literature 
we find a strong age dependence on how well SSP models (with underlying Kroupa or Salpeter-type stellar mass 
functions) fit the light-to-mass ratio of clusters. Based on this result we suggest that the large scatter in the 
light-to-mass ratio of the youngest clusters is not due to variations in the underlying stellar mass function, but 
instead to the rapidly changing internal dynamics of young clusters. Based on sampling statistics we argue that 
while W30 and G114 are extremely massive, they are consistent with being the most massive clusters formed in 
a continuous power-law cluster mass distribution. Finally, based on the positions of old globular clusters, young 
massive clusters (YMCs), ultra-compact dwarf galaxies (UCDs) and dwarf-globular transition objects (DGTOs) 
in K-space we conclude that 1) UCDs and DGTOs are consistent with the high mass end of star clusters and 2) 
YMCs occupy a much larger parameter space than old globular clusters, consistent with the idea of preferential 
disruption of star clusters. 
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1. Introduction 

Our concept of star clusters has changed rapidly dur- 
ing the past two decades. The first resolved young clus- 
ters with masses comparable to those of the traditional 
globular clusters (taken with the Hubble Space Telescope 
-Holtzman et al. 1992) confirmed the suggestions of 
Schweizer (1987) that mergers of galaxies may produce 
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'young' globular cluster sized objects. These results were 
rapidly followed by the discovery of additional young 
massive clusters (YMCs) in other galaxy mergers, as 
well as in dwarf, starburst, and normal galaxies (see re- 
views by Whitmore (2003) and Larsen (2004)). Even our 
own galaxy is producing YMCs with comparable masses 
and sizes to those observed in merging galaxies, e.g. 
Westerlund 1 (Clark et al. 2005). The apparent ubiquity 
of these objects has raised the question of how 'universal' 
their detailed properties are, in particular concerning their 
formation and subsequent evolution. 
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In order to address this and other questions regarding 
YMCs, we have begun a programme to obtain kinematic 
and structural properties of star clusters which lie at the 
extreme high-end of the distribution of observed (lumi- 
nous) masses. Our first result, presented in Maraston et 
al. (2004) was for the extremely luminous star cluster, W3, 
in the galactic merger remnant NGC 7252. Combining 
the velocity dispersion measured with UVES on the VLT 
(45 ± 5 km s"i) with the size determined from HST im- 
ages (Reft — 17.5±1.8) led to the dynamical mass estimate 
of 8(±2) X 1O''A40. This mass and that estimated from 
photometric methods (see Maraston et al. 2004 for details) 
were in excellent agreement, arguing that the stellar mass 
function within W3 was Salpeter-like. 

Using similar techniques as were employed in the above 
work, some studies have suggested that the stellar mass 
function in star clusters can vary substantially (e.g. Smith 
&; Gallagher 2001), while others have reported standard 
Kroupa (2002) or Salpeter (1955) type stellar mass func- 
tions (e.g. Larsen et al. 2004). These discrepant results 
have left the question of the variance of the stellar mass 
function of massive YMCs open to debate. 

Additionally, detailed knowledge of the internal dy- 
namics and structural parameters of YMCs has allowed 
a comparison between them and other gravitationally 
bound systems. In Maraston et al. (2004) we showed that 
W3 is too diffuse for its mass when compared with old 
globular clusters, whereas it is too compact relative to 
dwarf galaxies. However, we also showed that WS's prop- 
erties were extremely similar to those of massive point-like 
objects discovered in Fornax (Hilker et al. 1999). Based 
on this similarity, Bastian et al. (2005b) have suggested a 
mechanism which may allow massive star clusters to exist 
far from the main body of the host galaxy, namely the for- 
mation of massive clusters in the tidal debris of galactic 
interactions /mergers. 

Following up on these results, we have obtained high 
resolution UVES optical spectra of two additional highly 
luminous star cluster candidates in galactic merger rem- 
nants. The first is the cluster W30, the second bright- 
est cluster in NGC 7252. W30 has an estimated age of 
^ 300 — 500 Myr and a metallicity between half solar 
and solar, estimated from optical spectra (Schweizer & 
Seitzer 1998) as well as from optical and near-infrared 
photometry (Maraston et al. 2001). The observed magni- 
tude of W30 is my (W30)= 19.46 mag and it has a. (V - 1) 
colour of 0.63 mag (Miller et al. 1997). The second cluster 
in this study is G114 in NGC 1316. Based on optical and 
near-infrared photometry along with near-infrared spec- 
troscopy, Goudfrooij et al. (2001a,b) estimate an age of 
3.0 ± 0.5 Gyr for G114. This cluster, the brightest one in 
NGC 1316, has an observed magnitude mB(G114)= 19.63 
mag and {B — I) colour 1.87 mag. Throughout this work 
we adopt the distances to NGC 7252 and NGC 1316 
which were used in Maraston et al. (2004) and Goudfrooij 
et al. (2001a, b) respectively, namely 64.4 and 22.9 Mpc. 
Figures [H & |3 show the HST/WFPC2 planetary camera 
chip images of NGC 7252 and NGC 1316, respectively. 



along with the slit sizes and positions used. In Fig. ^ we 
also mark the massive star cluster W3. 

The luminosities of G114 and W30 ( assuming a 
Salpeter or Kroupa-type stellar IMF) imply that they have 
extremely high masses, more than 3 times that of the most 
massive globular cluster in the Galaxy, wCen which has a 
mass of - 3 X 10*^ Mq (Meylan & Mayor 1986). 

In this work we investigate the structural and kine- 
matic properties of these two massive star clusters in order 
to determine how they relate to "normal" young and old 
globular clusters. In §[5] we measure the effective radii of 
the two clusters using high-resolution HST imaging, and 
in § 121 we determine their velocity dispersions. We com- 
bine these results in § ^ to estimate the dynamical mass 
of the clusters. In §[3 we discuss the implications of our 
results in terms of the underlying stellar mass function 
of the clusters, and compare their properties with other 
bound stellar systems. Finally, we summarise the results 
and present our conclusions in § 

Fig. 1. i^555M^ image of the centre of NGC 7252 showing 
the slit size and position, as well as cluster W3 for refer- 
ence. The image is 740 pixels on a side which corresponds 
to a hnear distance of ~10.5 kpc at the assumed distance 
of 64.4 Mpc. 



Fig. 2. F450W image of the centre of NGC 1316 showing 
the slit position. The image is 800 pixels on a side which 
corresponds to a linear distance of '^4 kpc at the assumed 
distance of 22.9 Mpc. North is down and east is to the 
right in this image. 



2. Determination of the cluster sizes 

Structural parameters for clusters G114 and W30 were 
measured on Wide Field Planetary Camera 2 (WFPC2) 
and on Advanced Camera for Surveys (ACS) HST im- 
ages. The images for W30 are presented in detail in Miller 
et al. (1997). Here we note that W30 is located on the 
Planetary Camera chip. For G114 we measured the size 
on both the WFPC2 (data presented in Goudfrooij et 
al. 2001b) and ACS (presented in Goudfrooij et al. 2004) 
images. 

Sizes were found using the ISHAPE routine of Larsen 
(1999). This routine convolves the PSF with a specified 
model profile of varying sizes and fits it directly to the 
images. The outputs of this routine for the best fitting 
model are the FWHM of the major axis, the minor to 
major axis ratio, and the goodness of fit. For the WFPC2 
images we used a PSF generated by Tiny Tim (Krist & 
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Hook 1997) at the exact location of the cluster on the 
chip. We used the drizzled ACS images to eliminate the 
geometric distortion of the camera. The PSF for each filter 
of the ACS images was constructed using sources from 
ACS observations of the globular cluster 47-Tuc. 

For each cluster we fit multiple profiles to the images 
in each of the observed filters, including King, Moffat, and 
Gaussian profiles. 

2.1. NGC 1316:6114 

Using the FibOW WFPC2 observations, we fit eight pro- 
files to the image of the cluster. The results are shown in 
Tabled In addition to the set profiles, we also fit a Moffat 
profile with a variable index. The best fitting model is 
one with an index 1.76. Using this profile we measure an 
eflFective radius of 4.08 pc, which is remarkably close to 
the mean of the size determined using all the other pro- 
files (4.10 pc). We estimate the error on the size as the 
standard deviation of the size measured for the different 
profiles, 0.25 pc. 

We have also measured the size of G114 on Fb^bW 
and F814PF HST-ACS images. In Tabic d we show the 
results for the fits using the best fitting profile from the 
K^fPC^ observations. We find that the size determined on 
the F555Wacs images is ~ 15% larger than that found 
on the WFPC2 images. However the size measured on the 
FS14Wacs images is 15% smaller than on the WFPC2 
images. 

As the three images all give approximately the same 
result we conclude that G114 is resolved, and we assign 
the size of i?e//(G114)= 4.08 ± 0.55 pc. 



Table 1. Effective radius of NGC 1316:G114 for various 
models as determined from the F450W^ WFPC2 image. See 
text for details. The best fitting model (lowest y^) is shown 
in bold. 



WFPC2 


Type 


index 


-Rotf(pc) 






Gauss 




4.34 


1.93 




King 


5 


4.27 


1.54 




King 


15 


3.89 


1.13 




King 


30 


3.66 


1.24 




King 


100 


4.14 


1.87 




Moffat 


1.5 


4.34 


1.09 




Moffat 


2.5 


4.09 


1.14 




Moffat" 


1.76 


4.08 


1.0 


ACS 


Filter 


Type 


index 






F555W 


Moffat 


1.76 


4.83 ±0.25 




F814VF 


Moffat 


1.76 


3.68 ±0.1 



^ Best fitting model of the WFPC2 observations. 
^ Calculated individually for the different cameras. 



2.2. NGC 7252:W30 

Using the FbbbW WFPC2 image of W30, we find that 
the best fitting model is a King profile, with concentra- 
tion factor 100. Using this profile we measure an effective 
radius of 8.3 pc. However the mean size of the other pro- 
files gave 10.3 pc. Table|21gives the determined size of W30 
for different profiles. There is a clear decreasing trend 
as one goes to smaller radii (i.e. closer to the best fitting 
radius). However, a King profile with concentration factor 
30 along with a Moffat profile with index 1.5 are also ac- 
ceptable fits in terms of their x^. We therefore estimate 
the effective radius of W30 to be 9.3±1.7 pc, where the 
error of 1.7 pc convers the full range of acceptable profile 
fits. 

Table 2. Effective radius of NGC 7252 :W30 for various 
models as determined from the F555W^ WFPC2 image (PC 
chip). See text for details. The best fitting model (lowest 
X^) is shown in bold. 



WFPC2 Type 


index 


-Rcft(pe) 


2 / 2 


Gauss 




13.2 


3.67 


King 


5 


12.5 


2.93 


King 


15 


11.0 


1.82 


King 


30 


10.2 


1.31 


King 


100 


8.3 


1.00 


Moffat 


1.5 


11.0 


1.21 


Moffat 


2.5 


12.0 


2.19 



3. Velocity Dispersion 

We observed NGC 7252:W30 and NGC 1316:G114 with 
the Ultraviolet Echelle Spectrograph ( UVES) mounted 
on the ESO/VLT on the nights of Sept. 13-16th, 2004. 
We used the red arm, CD#3 grating centred on 5200A. 
This resulted in a wavelength coverage from 4200A to 
6200 A and a resolution of « 5 km s'^ at 5200 A. The 
data were reduced and extracted using the on-line UVES 
pipeline with the relevant bias and flat-frames. Cosmic- 
rays were also removed using the pipeline. After extrac- 
tion, each spectrum was corrected to the helio-centric ve- 
locity frame, and summed to create the total spectrum 
for each cluster. The total exposure times were 8.67 hours 
on G114 and 25.1 hours on W30. Figures □& H show 
the positions and lengths of the UVES slits for clusters 
W30 and G114 respectively, superimposed on 775 T" plan- 
etary camera images. We note that the background near 
the positions of each of these clusters is devoid of spuri- 
ous sources, which allowed a clear background subtraction 
(determined using a spline function) at the position of the 
clusters. 

Additionally, we observed several template stars 
which were used to complement our existing tem- 
plate catalogue (Maraston et al. 2004), namely 
HR 35 (F4 V), HR 8709 (A4 V), HD 203638 (KO HI), 
and HD 212574 (A6 V) where the designation in the 
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brackets refers to the spectral type of the star. The stars 
were reduced in the same way as described above. 

The determination of the velocity dispersion of each of 
the clusters was carried out in exactly the same way as was 
done for W3, which is described in detail in Maraston et 
al. (2004). In summary we used an adapted version of the 
Fourier Correlation Quotient (FCQ, Bender 1990) method 
as implemented by Bender et al. (1994), using templates 
chosen to match the stellar populations within each clus- 
ter. The templates were chosen to have temperatures and 
gravities (i.e. luminosity classes) appropriate to stars at 
the main-sequence turn-off point and giants stars in syn- 
thetic stellar populations of the same age as each clus- 
ter. Their contributions to the composite template were 
weighted using the weights predicted by the SSPs at the 
appropriate age (Maraston, 2005, e.g. her Fig. 13). 

3.1. NGC 1316:G114 

The results of the determination of the velocity dispersion 
for G114 is shown in TableEl Due to the high S/N ratio of 
the data and the large number of metal lines in the optical 
part of the spectrum (due to the dominance of cool stars at 
the cluster's age of 3 Gyr) we were able to use the full 
spectral range to determine the velocity dispersion. We 
note that the results do not depend crucially on the as- 
sumed stellar template. The adopted one-dimcnsional ve- 
locity dispersion for G114 is ct(G114) = 42.1±2.8 km s-\ 
which is the average over the full wavelength range of 
Templates 2 and 3, which should be the closest match 
to the actual stellar population based on the models of 
Maraston (2005). Fig. shows the blue section of the ob- 
served spectrum of G114 (black), the best fitting broad- 
ened stellar template (red) and the difference between the 
two (green) . All spectra shown in this work have been di- 
vided by the continuum and had a value of one subtracted 
from them, to place the average value at zero. 

We have also measured the heliocentric line-of-sight 
velocity of G114, 

w(G114) = 1292 ±3 km s'^ 

This is consistent with the results of Goudfrooij et 
al. (2001a), who measured 1306 ± 26 km s^^. 

3.2. NGC 7252: W30 

Based on comparisons between the optical/near-infrared 
colours (Miller et al. 1997, Maraston et al. 2001) and op- 
tical spectroscopy (Schweizer & Seitzer 1998), W30 and 
W3 appear to have very similar ages and metallicities. 
Because of this, we have used the same stellar template 
to determine the velocity dispersion of W30 as we used 
for W3 (Maraston et al. 2004). Similarly as was done 
for W3, we have limited our analysis to the region red- 
ward of H/?, in particular concentrating on the region 
around the Mg triplet (AA5167, 5172, 5183A) and Fe lines 
(AA5270,5335A). Table H shows the results of the fitting 



on specific regions of the spectrum. We adopt the value 
(t{W30) — 27.5 ±2.5 km s^^ which is the average between 
the value found for fitting solely on the Mg lines and fitting 
on the full region of interest (between 5150A and 5350A). 

Fig.0]shows the spectrum of cluster W30 in the fitting 
region (black), the best fitting broadened template (red), 
and the difference between the two (green). 

The measured line of sight velocity of W30 is 

v{W30) = 4614 ± 1 km s-\ 

which is in very good agreement with previous mea- 
surements, namely 4624 ± 17 km s^^ (Schweizer & 
Seitzer 1998). 

Table 4. Derived velocity dispersion for different portions 
of the spectra of W3G in NGC 7252. All values are given in 
km a~^. 





Fe" 


Full'= 


Adopted 


24.7 


30.3 


28.0 


27.5 ±2.5 



^ Just the Mg triplet, from 5156A to 5187A. 

^ Red- ward of the Mg triplet, mainly dominated by Fe lines, from 
5218A to 535lA. 

The full spectrum used, from 5155A to 5324A. 



4. Dynamical Masses 

Assuming that the clusters are in virial equilibrium, we 
can estimate their virial masses using the relation 

Mdyn = (1) 

(Spitzer 1987, p. 11-12) where 77 is a dimensionless param- 
eter which depends on the cluster profile adopted. Here 
we adopt 77 — 9.75, which is approximately valid for most 
of the globular clusters in the Milky Way. However, we 
will return to this adopted value, and its consequences in 
§ 15.1.11 Inserting the numbers derived above we find that 

A^d.vn(G114) = 1.64(±0.13) X lO^TW© 

A^dy„(H^30) = 1.59(±0.26) x lO^X© 

where the errors do not include uncertainties in the dis- 
tance to the clusters. The adopted properties of W3, W30, 
and G114 are shown in Table [3 

We note, however, that the value of the constant in 
Eq. n known as 77, may change as a function of age of 
the cluster (Boily et al. 2005). This effect will be most 
dramatic in the first 30 Myr of a cluster's lifetime, and 
the size of the effect depends on the surface mass density 
of the cluster in the sense that the clusters with the highest 
densities will be the most affected. We will return to this 
point in Section [5. 1.11 
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Fig. 3. The observed (observed divided by the continuum minus 1 and smoothed by 15 pixels) spectrum of G114 
is plotted in black. The red shows the best fitting template, broadened by 42.1 km s^^. The green is the difference 
between the two, shifted downward by a constant for clarity. This is just one region that we fitted, Table 13 shows 
all the regions that were used. The large scale undulations (~ 30A) are due to the response of the individual orders 
of the echelle spectra. Since the scale length of these variations is much larger than those of the spectral features of 
interest (< 5 A) they do not affect the derived velocity dispersion. 



Table 3. Velocity dispersion measurements of NGC 1316:G114. 



Template Spectral Type cr (km s ^) 

4360-5115A 5245-5432A 5486-5855A 5943-6113A Average 



HR 35 F4 V 41.8 47.3 43.9 43.4 44.1 
HD 203638 KO III 43.6 41.9 48.0 47.6 45.3 
Template 1" - 39.9 36.3 41.0 45.6 40.7 
Template 2'' - 40.3 37.8 40.7 44.7 40.9 
Template 3" - 41A 41^ 45^2 45S 43.4 



" 50% KO III star and 50% F4 V star 
^ 60% KO III star and 40% F4 V star 
" 70% KO III star and 30% F4 V star 



Fig. 4. The observed (smoothed by 5 pixels) spectrum of W30 is plotted in black. The red shows the best fitting 
template, broadened by 28.0 km s~^. The green in the difference between the two, shifted downward by a constant for 
clarity. The top two panels show the entire wavelength range which was used in the fitting, while the bottom figure 
shows a blow up of the region including the Mg triplet, which are the strongest lines in this region. 



These results confirm the extremely large masses of 
G114 and W30. However, it is important to put these 
clusters in the context of the mass functions of the full 
cluster systems of their respective galaxies in order to test 
to what extent they really are outliers. ^ 

To do this we have performed a series of monte-carlo 
tests of the mass function of cluster populations. We as- 
sume an initial power-law mass distribution of the clusters 
within each galaxy of the form N{M)dM oc M~°'dM, 
with a ~ 2 (e.g. Miller et al. 1997) and also that this dis- 
tribution gets filled randomly. In the case of the NGC 7252 
system, W3 (the most massive cluster in the NGC 7252) is 
~ 5 times more massive than W30 (the second most mas- 
sive cluster in this system). Under these conditions, we 
expect that ~ 20% of the realizations of the cluster pop- 
ulations will have a factor of five or greater between the 
most massive and the second most massive clusters within 
the system^. The third brightest cluster in NGC 7252 sys- 
tem (W6) is ~ 0.2 magnitudes fainter than W30, corre- 
sponding to a mass difference of only ~ 17% (assuming 
a common age). Hence, we see that while W3 and W30 
are extremely massive clusters, they are compatible with 
being the most massive clusters of a continuous power- 

^ As noted in Miller et al. (1997) and Goudfrooij et 
al. (2001a,b) both NGC 7252 and NGC 1316 have large pop- 
ulations of luminous star clusters which appear to be coeval 
with the merger of the host galaxies. 

^ We have also implicitly assumed that the ratio between the 
lowest mass cluster and the highest mass cluster is <^ 0.001. 



law distribution, which we note continues to the detection 
limit. 



The intermediate aged cluster population of NGC 1316 
can also be readily explained by the same argument as 
above, as the second and third most massive clusters in 
this system (G114 is the most massive) are only 0.47 and 
0.54 magnitudes fainter respectively. This corresponds to 
less than a factor of two in the luminosity (and mass 
assuming that the clusters have similar ages and metal- 
licities). A difference of two or greater in the ratio of 
the most massive and second most massive clusters was 
found in ^ 50% of the realizations of a cluster population. 
Goudfrooij et al. (2004) have reported that the bright end 
of the luminosity function (which we assume to represent 
the mass function) is well approximated by a power-law 
of the type used above. 



Since NGC7252:W3, NGC7252:W30, and 
NGC1316:G114 can be readily understood through 
sampling statistics we will assume that they are simply 
the most massive clusters of a continuous distribution. 
In § 15.21 we will use this assumption and the detailed 
properties of these clusters to understand more enigmatic 
objects, namely the dwarf galaxy transition objects 
(DGTOs) and ultra-compact dwarf galaxies (UCDs). 
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Table 5. The properties of the massive star clusters. 



Name 




Roff 




(T 




age 


Mass 




(mag) 


(pc) 


(km 






Gyr 




NGC 7252:W3 


-16.2 


17.5 ± 1.8 


45 


± 


5 


0.4 ± 0.15 


8 ± 2 


NGC 7252:W30 


-14.6 


9.3 ± 1.7 


27.5 


± 


2.5 


0.4 ± 0.15 


1.59 ± 0.26 


NGC 1316:G114 


-13.0 


4.08 ± 0.55 


42.1 


± 


2.8 


3± 1 


1.64 ± 0.13 



^ We have corrected for foreground extinction of NGC 7252 and NGC 1316 of Av — 0.03 and Ay — 0.0 respectively. 



5. Discussion 

5.1. Stellar mass functions 

A comparison between the mass of a cluster derived using 
photometric methods and the mass determined through 
kinematic arguments, allows an independent test of the 
assumptions that went into each estimate. The assump- 
tion that has garnered the most attention in recent years, 
is that of the underlying stellar mass function (MF) of 
the cluster. In order to estimate the photometric mass of a 
cluster, mass-to-light ratios from simple stellar population 
(SSP) models are required. These, in turn, are heavily de- 
pendent on the assumed stellar mass function. Therefore, 
assuming that all the other assumptions are valid (such as 
the state of equilibrium, correct extinction determination, 
and stellar evolutionary tracks) any difference between the 
mass of a cluster derived in these two ways is caused by 
the difference between the input stellar mass function and 
that of the cluster. 

Studies that have used this technique have re- 
ported strongly divergent results. For example. Smith & 
Gallagher (2001) and McCrady et al. (2005) have reported 
that the YMC M82F is deficient in low-mass stars, rela- 
tive to the standard Salpeter (1955) (a single power-law 
mass distribution from the lower to the upper mass limit) 
or Kroupa (2002) type MFs (a single power-law above 
1A40 and significantly flatter below this limit). However, 
Maraston et al. (2004), Larsen et al. (2004), and Larsen & 
Richtler (2004) have shown that clusters in a wide variety 
of galactic environments are consistent with a Salpeter or 
Kroupa- type MF. 

We therefore carry out this experiment for the two 
massive clusters W30 and G114. In Fig. we show the 
light-to-mass ratio from the Maraston (2005) models for 
solar metallicity and a Salpeter (dashed line) & Kroupa 
(solid line) stellar mass functions. Over-plotted in red are 
the two clusters in the present study as well as W3 from 
our previous study. These three clusters all lie impressively 
close to the value using a Kroupa MF. Hence they are most 
likely not deficient in low mass stars. 

5.1.1. Cluster measurements from the literature 

In order to compare our results to other young clusters, 
we have taken a sample of clusters with velocity dispersion 
and radii measurements from the literature. The clusters, 
their parameters, and the corresponding references are 



listed in Table We have taken the fundamental param- 
eters (age, extinction, brightness, velocity dispersion, ra- 
dius, and distance modulus) directly from the given refer- 
ence. In some cases the V^-band magnitude was not given, 
although we note that all clusters with ages greater than 
20 Myr have observed V-band magnitudes. These older 
clusters will constitute the main part of our comparison. 
In those cases where V-band magnitudes were not avail- 
able we transformed the given magnitude to the V^-band 
using the colours in the Maraston (2005) SSP models at 
the appropriate age (which assume a Salpeter IMF). 

We have estimated the mass of each cluster using Eq.^ 
and used their V^-band magnitudes (and given ages) to 
place them in Fig. (blue points). From this figure it is 
clear that the amount of deviation from standard stellar 
mass functions (Kroupa or Salpeter-type) is heavily age 
dependent, with the older clusters (with the exception of 
M82F'^) all consistent with a Kroupa or Salpeter-IMF and 
the youngest clusters showing a large amount of scatter. 
Below we outline three possible explanations for this. 

A first possibility for the age-dependent scatter in 
Fig. El is that 77 (the parameter in the numerator of Eq.^l 
changes as a function of time (Boily et al. 2005). This is 
caused by mass-segregation in young clusters and further 
internal dynamical evolution of the cluster. The variation 
of 7] is expected to also be heavily dependent on the sur- 
face density of the star cluster, with higher surface den- 
sities leading to larger variations of rj. Fig. El shows the 
mean surface density within the half-light radius (using 
the estimated virial mass of the cluster) of star clusters 
vs. their measured velocity dispersions. The small blue 
filled triangles represent old globular clusters in our galaxy 
(McLaughlin & van der Marel 2005), the small green filled 
circles are globular clusters in M 31 (a collation of data 
from McLaughlin & van der Marel in prep.), the small 
upside-down magenta triangles are old GCs in NGC 5128 
(Martini & Ho 2004), and the large red circles are young 
massive star clusters in a variety of galaxies (listed in 
Table EJ. AH of the YMCs in Fig. El have surface densi- 
ties above lO'^AIopc"^ and most are above lO^A^opc"^, 

^ M82-F, located in Fig. |Klat log (age) = 7.8, may be a de- 
viant point in the diagram due to uncertainties in its age and 
extinction. Further studies to pin down the exact values would 
be desirable. Additionally, we note that its elliptical shape, 
crowded environment, and peculiar radial velocity suggest that 
the cluster may have been gravitationally influenced by its sur- 
roundings, and hence may not be in virial equilibrium. 
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which is the regime where i] is expected to vary strongly 
(Boily et al. 2005). 

A second possibility is that many of the youngest clus- 
ters are not in dynamical equilibrium. This could be due to 
external gravitational effects (e.g. close passages to mas- 
sive GMCs) . As clusters are born in gas-rich environments 
this is a likely possibility. A lack of equilibrium could also 
be caused by the rapid expulsion of the gas left over from 
the star formation process (assuming a non-100% star for- 
mation efficiency). This can have a severe influence on 
a young cluster (e.g. Boily & Kroupa 2003). Bastian et 
al. (2005a) have suggested that rapid gas removal may 
be responsible for the dissolution of 70 — 90% of clusters 
within the first 10 Myr of their lives, independent of clus- 
ter mass. The lack of dynamical equilibrium is also sup- 
ported by the work of de Grijs, Wilkinson, & Tadhunter 
(2005) who showed that the clusters which deviate the 
most from the old globular cluster My vs. log(cr) relation 
are found in the highest density environments, and hence 
are the most likely to be affected by interactions with the 
external environment. 

Finally, a third possibility for the observed age- 
dependent scatter in Fig. [S] is that only clusters with 
Kroupa or Salpeter-type stellar mass functions survive for 
more than ^ 100 Myr. Smith & Gallagher (2001) suggest 
that if the cluster M82 F has a significantly top-heavy 
stellar IMF (i.e. truncated below 2-3 Mq) it will lack the 
gravitational potential to remain bound due to stellar evo- 
lutionary mass loss after 2-3 Gyr. However, in Fig. [S] the 
clusters which deviate the most from Salpter or Kroupa- 
type stellar IMFs have light-to-mass ratios below the ex- 
pected value. This implies that they are over-abundant in 
low mass stars relative to Salpeter or Kroupa-type IMFs. 
Since low mass, long-lived stars provide the gravitational 
potential to keep a cluster bound, we would expect these 
clusters to be long-lived and hence to see old clusters 
which also have light-to-mass ratios below that expected 
for standard IMFs. Such clusters are clearly lacking in 
Fig. |S1 arguing against this possibility. 

5.2. More clues from n-space 

5.2.1. The future evolution of W30 and G114 

Following the analysis by Maraston et al. (2004) we at- 
tempt to place W30 and G114 in the broader context of 
gravitationally bound stellar systems. For this we exploit 
the re-definition of the fundamental plane known as k- 
space (Bender et al. 1992) which combines the three fun- 
damental observable parameters (radius, surface bright- 
ness, and velocity dispersion) into physically motivated 
values. 

Figure |3 shows the position of W3 (asterisk), W30 
(upward triangle), and G114 (downward triangle) in the 
k1 - k2 plane. In this space k1 traces the mass of the 
system, while k2 measures the compactness of a system 
for a given mass (the product of M. /L and surface bright- 
ness). For comparison we also show the mean positions 
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Fig. 6. The measured velocity dispersion of star clusters 
versus their mean surface density within the half light 
radius. The blue triangles, small green circles, magenta 
upside-down triangles and the large red circles represent 
old globular clusters in the Galaxy, M31, NGC 5128, and 
young massive star clusters in a variety of galaxies, respec- 
tively. Note that the masses have been estimated assuming 
virial equilibrium. However, many of the young clusters 
have surface densities greater than 10^ A^o/pc^, which is 
the regime where rj is expected to vary significantly (Boily 
et al. 2005). 



of bulges and ellipticals (B-l-E), dwarf ellipticals (dE), 
M32 (all taken from Burstcin et al. 1997, and assum- 
ing Hq = 75 km/sec/Mpc). Next, we add the point- 
like objects in Fornax (the ultra-compact dwarf galax- 
ies or UCDs) (the average value of the four objects pre- 
sented in Drinkwater et al. 2003). We also plot the po- 
sitions of old globular clusters in the Milky Way (small 
blue triangles) (McLaughhn & van der Marel 2005), M31 
(small green dots) (a collation of data from McLaughlin 
& van der Marel in prep.), and NGC 5128 (magenta up- 
side triangles) (Martini & Ho 2004), for which we as- 
sumed a constant {B — V) colour for old metal poor GCs, 
namely 0.7 mag. We also add the young clusters taken 
from the literature (see Table EJ, which are shown as 
red points. The large filled green squares represent the 
Nuclear Star Clusters (NCs) in bulge-less disk galaxies 
(Boker et al. 2004, Walcher et al. 2005, see Table El for 
details). Finally, we add the dwarf galaxy transition ob- 
jects (DGTOs) found in the Virgo galaxy cluster (Ha§egan 
et al. 2005) as filled blue squares. For the DGTOs we as- 
sumed that {B—V)=l, typical of a 10 Gyr solar metallicity 
SSP. 

The arrows which begin at W3, W30, and G114 rep- 
resent the evolution of the cluster in this space when the 
clusters are 'aged' to a common age of 10 Gyr using the 
Maraston (2005) SSP models. Note, however, that the SSP 
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Fig. 5. The derived light-to-mass (l/-band) ratios as a function of the age of the clusters. Over-plotted are the 
ratios (for solar metallicity) from the Maraston (2005) SSP models with a Kroupa (solid line) and Salpeter (dashed 
line) stellar mass function. The clusters older than a few x 10^ Myr (except M82F) are well fit by the models, whereas 
some the youngest clusters show significant deviation. This may be caused by a strongly varying rj in the early stages 
of a cluster's lifetime (Boily et al. 2005). For M82F, we show the position of this cluster for two ages, 4 Myr and 
60 Myr connected with a dotted line. 



model tracks do not take mass loss (and hence fading) by 
evaporation or due to external perturbations into account. 

As was found for the massive cluster NGC 7252 :W3 
by Maraston et al. (2004), we find that NGC 7252 :W30 
and NGC 1316:G114 will evolve into the region of k- 
space occupied by the UCDs and DGTOs. This shows a 
strong similarity between the most massive star clusters 
and these enigmatic objects, and may suggest that they 
formed through similar mechanisms. 

Additionally, we can estimate the amount of mass loss 
which is expected to occur within W30 and G114. From 
the SSP models of Maraston (2005) we see that a clus- 
ter (assuming solar metallicity and a Kroupa stellar IMF) 
is expected to lose ~ 18% of its mass between the ages 
of 400 Myr and 10 Gyr (i.e. between the present age of 
W30 and the age of globular clusters). G114, with an age 
of ~ 3 Gyr is only expected to lose ~ 5% of its current 
mass to stellar evolution. Using the analytic expressions 
for mass loss in a tidal field of Lamers et al. (2005, eq. 6), 
we note that due to the strong dependence on cluster 
mass, neither of the clusters studied here are expected 
to lose a significant amount of mass due to disruption. 
W30 is expected to lose ~ 8% of its mass over the next 
10 Gyr while G114 is expected to lose just ~ 5% of its 
mass over the next 7 Gyr. For this calculation we have 
assumed (the average time for a lO^'M© cluster to dis- 
rupt) to be 1 Gyr, based on the galactic GC population 



(Boutloukos & Lamers 2003). However, we note that the 
conclusions reached are not significantly affected by the 
choice of 14. As these small changes would barely be visi- 
ble in Fig. Eland would only add confusion, we choose not 
to show them. 



5.2.2. The relation between YMCs and old globular 
clusters 

Figure [S] again shows the k1 — k2 projection of K-space, 
except with all of the YMCs (red points) aged to 10 Gyr 
using SSP models, again assuming only passive stellar 
evolution of the cluster. Here we see that many of the 
YMCs have evolved 'past' the globular cluster region and 
into the space occupied by W3, W30, the UCDs and the 
UGTOs. Burstein et al. (1997) suggest that the tightness 
of the GC relation in K-space may be due to the pref- 
erential destruction of star clusters outside a rather nar- 
row region of parameter space (e.g. mass and radius, see 
also Fall & Rees 1977). Along the same lines, Gnedin & 
Ostriker (1997) show that only a narrow region of the 
mass-radius plane of GCs is stable against disruption, and 
suggest that the initial parameter distribution may have 
been much larger than what is observed today. This may 
be what we are seeing in Fig. |H1 where the young clusters 
occupy a much larger region of K-space (in terms of mass, 
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radius, and compactness) than their older globular cluster 
counterparts. 

We note that many young star clusters are not ex- 
pected to survive for more than ~ 100 Myr (e.g. Bastian et 
al. 2005a; Fall, Chandar, & Whitmore 2005) due to inter- 
nal and external disruption mechanisms. Thus, it can be 
expected that many of the youngest YMCs in our sample 
will never survive to comparable ages as the Galactic GCs. 
We do not wish to imply in Fig.|Slthat all the YMCs will 
survive to comparable ages, only that once differences in 
their stellar populations are taken into account, the YMCs 
occupy a much larger parameter space than their old GC 
counterparts. In particular we note that YMCs tend to 
display extended envelopes in contrast to the tidally trun- 
cated older GCs (e.g. Schweizer 2004). The loss of such 
extended envelopes (i.e. due to interaction with their en- 
vironment) may significantly change the position of the 
young clusters in Fig. |H1 As NGC 1316:G114 is ~ 3 Gyr 
old it is likely that it has already lost its extended enve- 
lope, which may explain why it falls on the relation for 
old GCs in Fig. El 

The similarities between young massive clusters and 
old globular clusters have been shown in a number of re- 
cent works (e.g. Kissler-Patig 2004; de Grijs, Wilkinson, 
Tadhunter 2005). In particular, Kissler-Patig (2004) 
showed that YMCs will follow a very similar My-cr rela- 
tion (which is one projection of K-space) as old GCs once 
the fading of the young clusters (due to stellar evolution) 
is taken into account. However, the dangers of using such 
a projection can be seen when comparing Fig. 1 in Kissler- 
Patig (2004) with Figs.Q&IHlin the present work. In the 
My-cr projection, the UCDs appear to be quite consistent 
with the distribution of young and old clusters. However, 
in At-space (which also includes information on the size) we 
see that the UCDs are disjoint from the globular clusters. 
Additionally, as stated above, in K-space the young clus- 
ters clearly occupy a much larger parameter space than 
the old GCs, contrary to what is seen in the My-u projec- 
tion. Finally, we note that the UCDs, DGTOs, and YMCs 
with masses above a few x 10^ Mq all appear very similar 
with respect to their scaling relations (i.e. mass, veloc- 
ity dispersion, size, and mean mass density), suggesting 
a common formation mechanism (Kissler-Patig, Jordan, 
Bastian 2005). 

6. Conclusions 

We have presented velocity dispersion, effective radius, 
and hence dynamical mass measurements of two extremely 
massive clusters in galactic merger remnants. These re- 
sults confirm that galactic mergers can produce star clus- 
ters with masses well in excess of the most massive glob- 
ular clusters in the Milky Way. However, we have also 
shown that while these clusters are extremely massive, 
they are consistent with being the high-mass end of a con- 
tinuous power-law distribution of star clusters, suggesting 
that cluster formation is (mass) scale independent. 



Comparing the light-to-mass ratios of W30 (in 
NGC 7252) and G114 (in NGC 1316) to those predicted 
by simple stellar population models (at the ages of the 
clusters), shows that both of these clusters are consistent 
with having Kroupa-type stellar mass functions. Applying 
the same analysis to other young clusters taken from the 
literature shows that there is a significant age dependence 
on how well SSP models fit the light-to-mass ratios of 
young clusters. Therefore, it is possible that the deviation 
from the light-to-mass ratio of young clusters from that 
predicted by SSP models is not due to a varying stellar 
mass function, but instead reflects the state of equilibrium 
of the youngest clusters. 

We have shown that W30 and G114 currently reside 
at the high-mass tip of the old globular cluster distribu- 
tion in K-space (a re-definition of the fundamental plane) . 
Both clusters, along with many young clusters taken from 
the literature, are likely to evolve into the space occu- 
pied by the so-called ultra-compact dwarf galaxies (UCDs) 
and the dwarf galaxy transition objects (DGTOs). This 
shows that young massive clusters and UCDs/DGTOs 
share many similar properties and suggests that the enig- 
matic UCDs/DGTOs may have formed in a similar man- 
ner as the most massive globular clusters in mergers, i.e., 
under rather violent circumstances. 

Additionally, we showed that young massive star clus- 
ters will occupy a much larger region of K-space than 
presently occupied by old globular clusters. This is con- 
sistent with the interpretation that star clusters are born 
with a larger range of parameters (e.g. radius, mass, and 
compactness) than displayed by globular clusters, and de- 
structive processes whittle away at the initial full distri- 
bution with only clusters which have parameters within a 
small range surviving to old ages. 
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Table 6. The properties of young massive clusters taken from the literature. 
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For Fig. □ we have assumed (B — /) — 1.0 which corresponds to a simple stellar population of ~ 400 Myr and solar 
metallicity, although we note that the assumed {B — I) colour does not change the conclusions. 
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